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Abstract: 

Possible CP violating effects in Z — > 3 jet decays are investigated. The analysis 
assumes the presence of CP violating ZqqG couplings. The contribution of these 
couplings to the Z — ► qqG decay width is calculated for different cut algorithms, 
including nonzero quark masses. Various CP-odd observables are discussed and it 
is shown that their sensitivity can change significantly if one uses normalized or 
unnormalized momentum vectors for their construction. Optimal observables are 
proposed which allow to measure the new couplings simultaneously. 



1 Introduction 



The large number of Z bosons collected at the LEP collider allows detailed studies 
of Z decays. An interesting possibility is the search for CP violation beyond the 
one present in the Standard Model (SM). For general discussions and lists of further 
work we refer to |], p|. In this note we want to present new calculations for the Z 
decay rate which have already been applied in a recent experimental analysis |3|]. 
We further reinvestigate certain CP-odd observables and propose the measurement 
of optimal observables. 

The theoretical framework is defined in ref. p]. There, possible CP violating 
effects are parametrized by an effective Lagrangian, i.e. CP violating formfactors are 
assumed to be real and constant. This effective Lagrangian is being added to the 
SM and calculations are performed at tree level. Including operators of dimension 
6 and demanding SU(3)c x U(l) em invariance the effective Lagrangian relevant for 
Z — > 3 jet decays contains weak and chromoelectric dipole moments and a ZqqG 
vertex with couplings hy q and fiA q - 

The couplings hy q /A q correspond to dimension 6 operators which conserve the 
quark chirality. A possible way to generate such couplings is through loop diagrams 
involving non-standard Higgs scalars |3[], leading naturally to larger couplings for 
heavier quarks. In an analysis of Z decays it is therefore legitimate (albeit not 
compelling) to assume that the b quark couplings h Vb /Ab are much larger than the 
corresponding couplings for light quarks. In ref. such a scenario was used to 
explain the gap between the experimental and theoretical values of Rb with the 
existence of large CP violating couplings for the b. However, new (preliminary) 
data || seem to indicate that this discrepancy is no longer that significant. 

First experimental studies || have concentrated on the measurement of angular 
correlations of the jets in Z — > 3 jet decays. In the present framework these correla- 
tions are only sensitive to the couplings h Vq /A q (the chromoelectric dipole moment 
does not contribute and effects due to the weak dipole moment are suppressed by 
a relative factor m q /Mz). In order to obtain limits on the couplings hv q /A q from 
an independent measurement, the authors of used the differential two-jet rate, 
which requires the knowledge of the contribution of the new couplings to the decay 
rate. To avoid statistical uncertainties from event generators, analytic expressions 
for the rate are useful and shall therefore be presented in the following section. 

2 Anomalous contributions to the rate 

We consider the decays Z —>■ 3 jets on the Z resonance. At the parton level this is 
equivalent to the reaction 

e + e - ^ Z (p) - + q(k+) + G(k G ) , (1) 
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and we assume that the quark flavor (and perhaps also the jet charge) can be 
measured. It is convenient to describe this process in terms of production and decay 
matrices, p and R |jj . Taking the Z propagator in Breit-Wigner form and using the 
narrow width approximation, the total cross section factorizes, 

a(e + e- -> Z -> qqG) = a(e + e - -> Z)—Y{Z -> ggG) , (2) 

with the total Z width Y z and the on-shell condition p 2 = M§. The partial width 
r(Z — > qqG) is given as integral over the trace term of the decay matrix R. 
The phase space for the decay Z — > qqG can be written as 

]Iz -dy + dy_Q(y + y_(l - y+ - y_) - r 2 (y + + y_) 2 )6(cuts) . (3) 
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The (positive) integration variables y± are defined in terms of the q/q energies E± 
in the Z rest system as y± = 1 — 2E±/Mz; the gluon energy is Eg = Mz(y+ + y-)/2- 
We further defined r = m q /Mz- The second O-function implements the presence of 
cuts with cut parameter y, 

e(cuts) = 0(?/ + - - y)0(y +G - j/)6fo_ G - y) , (4) 

with 

f 2M^ 2 J B iJ B j (l - cos^) for JADE cuts, 
^ ~ { 2M^ 2 min{£ 2 , E 2 }(1 - cos for DURHAM cuts. 

Here i?^ is the angle between the momentum directions kj and kj = +, — , G) 
measured in the Z rest system. Note that the JADE algorithm agrees with the so- 
called 'E— cuts' yij = (ki + kj) 2 /M z only for massless fermions. These prescriptions 
can of course be applied to the decay Z — > ffj as well. 

In Fig. la) we show the shape of the phase space for massless fermions (the 
dashed triangle) and for b quarks (i.e. r = 1/20), both without cuts. Fig. lb) shows 
the effect of JADE and DURHAM cuts on the phase space for b quarks. In this case 
the border lines are nontrivial functions given implicitly by the zeroes of G(cuts) 
(although they look like straight lines). 

The SM result for the rate T(Z — > qqG) has been calculated to order a 2 in 0. 
Since the new couplings hy q /Aq are CP violating, their contribution to the width 
adds incoherently 0, 

T(Z - qqG) = + AT q , G . (6) 
Note that in the differential two-jet rate 

= T(Z Z qqG) Ty (7) 
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the anomalous part is not additive, since also the normalization is affected. The 
additional contribution to the decay width can be written as 



Al^-a = ^T up ( f + (y, r 2 )^ + h%] + /_(y, r 2 )[h 2 Vq - h\ q \ ) (8) 

with T uP = aM z /(24: sin 2 6 W cos 2 6w). The coefficient functions f±(y,r 2 ) are ob- 
tained from integrating the relevant parts of the trace term a of the decay matrix 
R (cf. eq. (|A.4j) ) over the phase space volumes shown in Fig. lb). Note that we 
consider only the contributions from the couplings hyq/Aq, i- e - we assume the dipole 
moments to vanish. Their contributions to the rate can be computed in a similar 
way with the full trace term a. 

The couplings h Vq /A q correspond to a pointlike interaction which is not present 
in the SM Lagrangian. Since there are no internal propagators, infrared singularities 
are absent and the anomalous contribution AT qq -G is therefore finite even if no cuts 
are applied. In this case one finds 
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/ + (0,r 2 ) = -(l r 2 - 29r 4 + 5r 6 + 30r 8 J y/l - 4r 2 



+ 48r 4 (l-r 2 + r 6 )ln 



1 + VI -4r 2 ' 



2r J 1 

2 

/_(0,r 2 ) = T — (3 + 26r 2 - 62r 4 + 60r 6 ) Vl - 4r 2 

1 + Vl -4r 2 ^ 



(9) 



+ 16r 4 (-l + 3r 2 -6r 4 + 5r 6 )ln^- 
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In the presence of cuts the coefficient functions f± can be computed analytically 
only for massless fermions. For nonzero fermion masses one can at least give a series 
expansion 

f±(y,r 2 ) = f£\y)+r 2 f£\y) + 0(r 4 ,rHnr 2 ) , (10) 

where the neglected terms for the heaviest fermion in question - the b quark - are of 
the order 4 ■ 1CT 5 . In ([10|) , f± (y) represents the exact result for massless fermions. 
For JADE cuts, the coefficient functions are 

fi°\y) = i( 1 - 3 ^ 2 -^ 1 - 3 ^ 4 -T^( 1 - 3 ^) 5 ' 

/f(y) = (1 - 3y) (-9 - + 13^-3^+8^(5-^) In ^ 1_y ^ 



2 v J (id 



f m (y) = 0, 

f ( -\y) = ^(l-s^ + ^i-sy) 4 
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For DURHAM cuts we obtain 



3y + Jy(8 + y) 



/| 0) ( y ) = —^24 — 350?/ - 720y 2 + 2570y 3 + 525y 4 - 2997?/ 

+ ^(400 - 754y - 321y 2 + 999^)^(8 + y) 

+ 8y 2 (-3 + 5?/ + 4?/ 2 -6?/ 3 )ln 

f { +\y) = \{ - 54 - my +ll(% 2 + 29?/ 3 + 351?/ 

+ 1 (120 - 380?/ + 5y 2 - 117y 3 ) yjy(8 + y) 

+ 8j/(- 13 + 20y - 3y 2 + 7y 3 ) In f 



3^ + ^(8 + 2/) 



(12) 



= o, 

/-%) = ^(6 + 28 2 /+132 Z / 2 - 322 2 / 3 + 35V 



|(- 68 + 122y-117y 2 )V2/(8 + y) 



+ 8y(-2 + 5?/-10y 2 + 7y 3 )ln f- 



3|/ + x/2/(8 + y) 



As a nontrivial check on these results we have /j.°' 1 ' ) (l/3) = since for the cut 
parameter y = 1/3 the phase space has shrunk to a point. 

In Fig. 2 we show the coefficient functions f±{y,r 2 ) of (|8|) for b quarks, i.e. for 
r = 1/20, both for the JADE and the DURHAM algorithm. For massless quarks, 
/ + (0, 0) = 2/5 whereas f-(y, 0) vanishes identically. For nonzero quark masses the 
first expansion coefficients in ( ]T0| ) are larger than the leading terms by a factor of 
~20. As a consequence, the mass corrections for b quarks amount to 5% (and not - 
as expected - to r 2 ~ 0.25%). 

For practical purposes one can neglect the contribution proportional to [h Vq — h 2 Aq ] 
in (|5p. Assume first that /_ vanishes. An experimental limit on the size of Ar 
translates then into a bound on (h Vq + h^) 1 ^ 2 which means that allowed couplings 

lie within a circle with a certain radius R in the hyq-hAq-p^snae. If /_ is turned on, 
this circle is deformed to an ellipse with half axes R/y/l ± e with e = /_//+. A 
numerical study shows however that 0.004 < e < 0.006 for all values of the cut 
parameter y, both for JADE and DURHAM cuts. 
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3 CP— odd observables 



The best way to measure CP violating couplings is of course with CP-odd observ- 
ables. Due to fragmentation effects, the parton spins cannot be reconstructed, and 
therefore CP-odd observables have to be built from parton momentum directions 
(the parton energies are determined once the angles are known [Q). This is in con- 
trast to tt production (see e.g. [§]) where the top spin can be traced through the 
decay. Observables that use only momentum directions, i.e. pure angular correla- 
tions, were proposed and investigated in [[[[]: 

= It^-^tI • C3) 




j 

where i, j denote Cartesian indices in the Z rest system. In addition we consider 
the following set of CP-odd observables: 

vf = (k + xL)«, (15) 

vf } = |k + xk_|(k + xk_) i , (16) 

if = (k + -L) 1 (k + xL) J + (mj), (17) 

if = |k + xk_|(k + -k_) i (k + xk_) i + (i~.7)- (is) 

These differ from fll3|), (|14T) only by factors of |k+ x k_| = sin-$ + _. In the effec- 
tive Lagrangian approach, the anomalous couplings hv q /A q correspond to pointlike 
interactions, which typically lead to event topologies which are not collinear. An 

additional weight factor sint^ suppresses collinear events and should therefore 

increase the sensitivity to hv q /A q - 

The observables listed above have definite transformation behaviour (as vec- 
tors and tensors), thus their expectation values are proportional to the Z vec- 
tor and tensor polarization s and s^. For unpolarized beams and taking the e + 
beam as 3-direction we have s = (0,0, 7 e i) and = |diag(— 1, — 1, 2) (here 7 e i = 
^QveQAe/ (gve + 9Ae) with the weak vector and axial vector Ze + e~ couplings gy e and 
gAe)- Thus what one really should measure are the components V3 and T 3 g . Also 
note that the tensor observables if do not change sign upon charge misidentifica- 
tion k + <-> k_, whereas for the vector observables the correct assignment of the jet 
charge is required, which makes them less useful experimentally. 

From now on we assume that only the b quarks have nonvanishing hy q /A q cou- 
plings. The expectation values of the observables should then be evaluated in an 
event sample containing Z — ► bbG decays whereas Z decays into light quarks can be 
used to test the CP blindness of the detector. 
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For the number N cut of Z — > events within cuts needed to see a 1 s.d. effect 
we have (writing generically O for one of the observables) 

(O 2 ) 

N ^ = joy* • ( 19 ) 

Following ref. JJ] we scale this number 

N= f^HG) N " (20) 

in order to give an estimate of the total number N of Z bosons (note however that 
this procedure disregards detector resolutions, b tagging efficiencies etc.). 

For the analysis of CP-odd observables we neglect terms quadratic in CP vio- 
lating couplings. The variances (O 2 ) are then determined from the SM amplitude, 
whereas the expectation values are linear functions of the anomalous couplings. In 
this approximation the sensitivity of the observable O defines a band in the \/N hyb~ 



N fiAb plane. In the case of the tensor observables the slope of this band is 
given by the ratio gvbl 9Ab of weak Zbb couplings, since the expectation values (T^) 
are proportional to the combination hb = hybQAb — ^Ab9vb [@] (the latter remains 
true also if nonzero dipole moments are allowed, cf. eq. ( |A.2| )). 

In Fig. 3 we show the result for the observables (|i3|)-([T8|). Plotted is the half 
width of the band in the y/~N hyb^y/N h A b plane as function of the cut parameter y. 
If y tends to zero, this width diverges for the observables V3 and T33 (long dashed 
curves). The best observables are V3 and T33 (solid lines) whereas the sensitivity of 



v} 3 ^ and T33 is slightly reduced. The gain in sensitivity when going from V3 ij (T< 
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to V3 (T33 ) can be easily understood: The unit vector n = (k + xk_)/|k + xk_| used 
in the construction of the observables V3 ^ and T33 gives only the direction normal 
to the decay plane, whereas the unnormalized vector k + x k_ has an additional 
information content, namely the angle between k + and k_. 

As a byproduct we have also recalculated the results of ref. for nonzero b 
quark mass (in mj = was used throughout) and found that the numbers iV 
given there increase typically by 3% - 8%. 



4 Optimal observables 

In this section we construct a set of optimal observables to measure the couplings 
hvb/Ab separately, following the procedure developed in || |10|| . Central assumptions 
are that one can neglect higher orders in the anomalous couplings and that the usual 
Gaussian error analysis applies. In linear approximation one writes the differential 
cross section for the reaction (|1]) as 

d a (e + e- ^Z^ bbG) = dcj)' (s + S v h Vb + S A h Ab + 0(h 2 Vb/Ab )) , (21) 
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where So represents the SM contribution. The phase space measure dtp' includes 
in addition to d<p of eq. (||) two integrations over angles relative to the e + beam 
direction p + . According to jlOfl , the optimal observables are then given by the 
ratios 

Sy/A 



O 



V/A 



S(\ 
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and any nonzero expectation value (Ov/a) 7^ signals CP violation. From the 
explicit formulae one easily derives 



Sn 



,(0) 

1 



</-■' +7d (p+k+rf } + (p+k_)«4' 



(o) 



1 1 




l + 2 


"(P + k + ) 2 - [ 



4 0) + 



+- 



4 0) + 



(p+k+)(p+k 



S 



V/A 



±f3 2 C F N c — f p + (k+x k_)\g Ab/Vb 



Ml 



+9Vb/Ablel 



y+ 
y- 



y- 
y+ 



1 ±4r 2 




(23) 



The term So contains only SM quantities; the relevant contributions a^°\ . . . , cjj to 
the coefficients of the decay matrix R can be found in The global factors j3, Cf, 
and N c are defined in the Appendix. 

The general formula given in [|10j for the sensitivity of these observables simplifies 
for the case of CP violating couplings. Since there are no contributions to the rate 



linear in h 



Vb/Ab, 



we have immediately Jd(f)' Sy = Jd(f)' Sa = 0. Moreover, only the 



trace term a contributes to the rate, i.e. / d<f>' S 
determined by the matrix 



Mi 



dcf) a 



-i 



= / d(j)' a. The sensitivity is then 



So 



(24) 



which defines - after scaling according to eq. (|2~0D - an ellipse in the \fN hy b - 
\fN liAb plane. In Fig. 4 we show the results of our numerical evaluation. The 
strong correlation in direction of the line /if, = is not surprising, since the other 
combination of couplings, hvbgvb — hAbgAb, is suppressed by a factor of j e \ ~ 0.16. 

(2,) 

In the same plot we compare the results for the best vector observable V% (long 
dashed) and the best tensor observable (short dashed). The optimal observables 
have a higher sensitivity, although the difference is only small. Note however that 
the construction of the optimal observables requires a measurement of the jet charge 
(cf. eq. (|2"3"D), similar to the vector observables. 
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5 Conclusions 



In this paper we have presented various calculations which can be applied in a search 
for CP violation in Z — » 3 jet decays. The effect of CP violating ZqqG couplings 
on the partial width T(Z —>■ qqG) was computed for different cut algorithms and 
for nonzero quark masses. An important result of our investigation of CP-odd 
observables was that the sensitivity of the conventional vector and tensor observables 
changes significantly depending on whether normalized or unnormalized momentum 
vectors are used for their construction. The analysis of optimal observables showed 
that the gain in sensitivity is only small. Nevertheless their experimental exploration 
might be useful in order to perform a complete analysis of the data. 
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Appendix 

The decay matrix R for the process Z — ■> qqG is defined as 

R ij = £<Z(eO|rt|<raG) (qqG\T\Z(ej)) , (A.l) 

where the summation extends over the q and q spins and the gluon polarisation. 
We work in the Z rest system, where Z(ei) denotes a Z with polarization vector 
along the i th Cartesian coordinate. The decomposition of R with the help of the q 
and q momentum unit vectors k_ and k + is given in 0. Note however that we did 
not normalize R with T^g and therefore we have | jd(f)Tr(R) = T^g (instead of 
|/#Tr( J R) = lini). 

The coefficients of that decomposition are calculated at tree level using in ad- 
dition to the SM Lagrangian the effective CP violating Lagrangian of |IJ which 

contains weak and chromoelectric dipole moments d q and d' q and the ZqqG cou- 
plings hy q and h^q- Contributions from the SM and terms linear in CP violating 
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couplings have already been given in [p]] . In order to make this article self-contained 
and to illustrate the different notations, we include the latter here as well: 



b [ s\y + ,y-) 



w(y+,y 



0. 



y+y- 



(k Vq g Vq -h Aq g Aq )(y + (l - y+) + y_(l -y-fj + 



+ i(rd q g Aq + r 2 (h Vq g Vq + h Aq g Aq )^j (y + + y_ 



(A.2) 



c [ i (y+,y-) = c { 2 ] {y + ,y-) = c 3 lj (y+,y-) 
c < i\y+,y-) = -c ( i\y-,y+) 







w(y+,y~] 



V- 



(I - y + ) 2 - Ar 2 . 



where 



w{y+,y-) 



hv q 9Aq — h Aq gvq , 

^Jy+y-i 1 -y+- y-) - r2 (y+ + y-Y 



| k + x k 



(A.3) 



Here gy q = T 3 — 2Q q sin 2 6 W and g Aq = T 3 are the SM Zqq couplings; for the 
definition of y± and r see eq. ([|). Out of the contributions quadratic in CP violating 
couplings we list only the trace term a^(y+, y~) and the vanishing coefficients: 



a {2) (y + ,y-) 



d r , 



(y+ + y- 



2y+y- 

} (y + y-) 2 
r iy+ + y-f 
2 y+y- 



2(y + + y-) + (l-4r 2 ) 



2, l ~y+-y- 



+ dqd'gvq 



3 - 2y+ - 2y- - r 



y+y- 
(y+ + y-f 



y+y- 



+ d q h Aq 
+ d' q g\ 

2 , o o % 1 



+ 2r(l +y + - 2r 



s(y+ + y-) 2 



y+y- 



^y+y- + r (~ 3 + y+ + y ) + r ' 



{y+ + y-f 
y+y- 



(A.4) 



+ d'a (9v q + 9Aq)~(. 1 ~ 2/+ - 2/-) + d'h Aq g Vq r(y + + y. 



+ d' q h Vq g Aq 



r r 3 (y + + y_) 3 



+ (^,-^,)j(3/h 



+ (h 2 Vq +h 2 Aq )^(4y + y^ + (y 2 + y 2 _ )(1 - y + - y_)) 
4 2) (l/+> 1/-) = 4 2) (j/+>£/-) = 4 2) (y + ,?/_) 



. 
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The nonvanishing coefficients b^\ &| , c[ , c^, and c 3 which do not contribute to 
the rate but are relevant for other (CP-even) observables have been calculated as 
well and are available on request. 

All coefficients in (|A.2|) , (|A.4j) have to be multiplied with the global factor 
(3 2 CfN c , with the coupling (3 = eg s / (sin 9w cos 9w) and the color factor Cp • N c = 
I • 3. The coefficients can of course be applied to the process Z — > f as well. 
In this case one has to replace the chromoelectric dipole moment d' q with the elec- 
tric dipole moment —df and the ZqqG couplings h Vq /A q with the Zffj couplings 
—fvf/Af (the minus signs stem from the normalization in [[[[). Moreover, the global 
factor is obtained from replacing in (3 the strong coupling constant g s with eQf and 
changing the color factor Cf-N c to 1-1 (1-3) for leptons (quarks). 
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Figure Captions 



Figure 1: Allowed phase space regions (Dalitz plots) for the decays Z — > ffV 
(V = photon or gluon) in the y+-y~ plane, a) Without cuts, for massless 
fermions (dashed) and for b quarks (solid), b) With cut y = 0.1, for b quarks 
with JADE (solid) and DURHAM algorithm (dashed). 

Figure 2: The coefficients f+(y, r 2 ) and f-(y, r 2 ) in the decomposition (^) shown 
as function of the cut parameter y for JADE (solid) and DURHAM cuts 
(dashed). The parameter r = m q /M z is put to 1/20 (6 quark). 

Figure 3: a) The half width of the sensitivity band in the \/N hvb~VN hAb plane 
as function of the cut parameter y. Shown are the results for the vector 
observables Vg (long dashed), V3 (solid) and Vg (short dashed), b) The 
same for the tensor observables T33 (long dashed), T33 (solid) and T33 (short 
dashed). 

Figure 4: Sensitivity contours of different observables in the \/N hyb^VN h,Ab 

plane. The results for the best vector observable (long dashed) and the 

(2) 

best tensor observable T 33 (short dashed) are compared to the ellipse from 
the optimal observables Oy, Oa (solid). 
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Figure 1 
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